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Abstract 



This article examines the decoherence of a macroscopic body us- 
ing a simple model of the environment and following the evolution 
of the pure state for the whole system. We found that decoher- 
ence occurs for very general initial conditions and were able to 
confirm a number of widely accepted features of the process. 
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1. INTRODUCTION 



There is growing interest in decoherence, particularly as it promises to 
deepen our understanding of the macroscopic world in terms of quantum 
principles [Q . If we suppose that S is a macroscopic body (say a point parti- 
cle with a very large mass M and position co-ordinate at the centre of mass) 
interacting with an environment E, we know from a host of examples that 
if we start with a factorized state of the two systems, the environment state 
vectors associated with the body's distinct positions (at least macroscop- 
ically) rapidly become orthogonal to each other, irrespective of the initial 
conditions. The spatial correlations of B, whose reduced density matrix be- 
comes diagonal on a position basis, are suppressed by decoherence. Thanks 
to this apparent collapse due to interaction with the environment, we see the 
classic characteristics of the macroscopic world emerge. During the above 
mentioned process temporal evolution also affects body B. In a number of 
examples, however (some in a seminal paper by Joos and Zeh[0]), wave packet 
spreading is ignored. This is justified both by the speed of the decoherence 
process and because mass M is large. We will adopt this approximation by 
formally ignoring the term of the body's kinetic energy in the total Hamilto- 
nian. This is rigorous with the limiting condition that M is infinite. Realistic 
sources of decoherence are described in the literature on the subject, such 
as scatteringl^, ^, Q and quantum gravity[^, Q but the focus is generally on 
useful and practical models such as the so-called Caldeira- Legett environ- 
ment 1^, which in any case reflects certain physical situations. Omnes 
recently put forward a general theory of the decoherence effect encompassing 
both Caldeira-Legett's harmonic model and the external environment con- 
sidered by Joos and ZehjQ, as well as being related to the quantum state 
diffusion model|§]. In the applications mentioned, the usual approach is to 
study the temporal evolution of the reduced density matrix we are interested 
in, which obeys an irreversible master equation[0]. Resorting to this tech- 
nique is almost inevitable as the overall system has many (or infinite) degrees 
of freedom. However, in my opinion at least, it would be useful to propose 
a further example of decoherence within the framework of an environment 
model using Schroedinger evolution of the initial pure state. Since the sem- 
inal paper by H. D. Zeh, published in the very first issue of Foundations 
of Physics|ll]], our understanding of the decoherence process has gradually 



2 



shifted from qualitative to quantitative. Our aim is to use pure state evolu- 
tion to obtain analytical results independent of a particular master equation. 
While there will be no ground breaking results, a new example could provide 
a further confirmation of our convictions, as well as offering a number of 
points worthy of consideration. This means the paper is didactic to a certain 
degree, but also represents a small step forward in our understanding of the 
decoherence process in quantitative terms. 



2. THE MODEL FOR THE ENVIRONMENT 

The example we will discuss concerns a one dimensional system in which 
the external environment is represented by n identical particles with x axis 
CO- ordinates Xi, 0:2, . . . , a;„. The macroscopic body's position is indicated by 
X. The Hamiltonian for the environment is based on the Hepp- Coleman [p!2| 



model, which appeared some time ago and was also discussed by Bell[0 
It was used to discuss the problem of measurement in terms of apparent 
collapse. (It was actually adopted to describe the system to observe, whereas 
here it will be used to describe the environment which acts as observer). The 
model considers the particle's kinetic energy as proportional to its momentum 
rather than as given by the customary quadratic term. The wave packet of 
a free particle of this sort moves along the axis in a single direction without 
changing form. Our Hamiltonian for the environment is 

He = a^pi + V{Xi,X2,. ■ ■ ,Xn) (l) 

where pi is the momentum, a is a constant quantity with the dimension of 
a velocity, and the potential term V{xi,X2, ■ ■ ■) includes both one particle 
interactions and two or many particle interactions. The expression assumed 
for the kinetic term is rather unconventional and the hamiltonian is consid- 
ered to be an appropriately idealized toy model. However, the intuition can 
be supported by reference to physical models. In the case of a free particle, 
for example, it may be considered as the caricature of an electron in one 
dimensional motion. Let us consider the ultra-relativistic approximation fre- 
quently used for very energetic electrons, and which amount to take the rest 
mass to be zero in Dirac equation. Then this transform in a couple of 
Weyl equations given by the hamiltonians = ±c cr • p applied to a two 
component spinor. As usual cr denotes the Pauli matrix in vectorial form. 
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Considering if_ and assuming that the spin of each particle in the one di- 
mensional system is antiparallel to the x axis (i.e. by considering only the 
component with Sx = —1/2 for all particle spinors) we obtain the kinetic 
part of He, with a = c. The free waves with positive momentum represent 
a massless particle with spin pointing in the opposite direction (a neutrino). 
The free waves with negative momentum, on the other hand, have negative 
energy and opposite helicity and should be referred to as antiparticles. The 
problem now is the other term in our environment hamiltonian, as it is not 
easy to imagine an interaction depending on positions, which in any case 
should be given in a relativistically correct form for particles of this type. 
We may also consider the electrons in the periodic field of crystal atoms, 
which may be thought of as free "dressed electrons". We know|T5| that for 



an infinite simple cubic lattice the one dimensional motion energy curves of 
these free "quasi particles" can be described with a good degree of accuracy 
in the form: 

W{q)n = An- BnCos{aq) 

where a is the lattice spacing and q the reduced wave vector. Interest is 
usually focused on the neighbours at the top and bottom of the band. But 
for our purposes we take one band (for example the first) and develop the 
dispersion relation above near qa = tt/2, assuming the zero of energy scale 
to be Aq: 

Wiq)o -Ao = Boiaq - n/2) = B^aq' = E{q') 



with the limitation \aq'\ < 1. We then follow Wannier's method [|I6| and 
replace q' with —id/dx in the last expression to obtain the effective hamilto- 
nian: 

H = E{-id/dx) = iBoa/h)p . 

Our model can therefore be thought of as describing a one dimensional system 
of interacting quasi particles in a limited pseudo-momentum zone. 

For the interaction of the environment with body B we will assume the 
linear coupling W^'^^^ = kX^Xi, where /c is a constant, already considered 
in the literature. It is similar, for example, to von Neumann's measurement 



interaction [0. This coupling is debatable in that it is not invariant under 
translation ||18|. For our purposes we will use it in this form, but will then 
show that it can safely be replaced by a harmonic coupling. If we assume B is 
free (apart from interaction with the environment) and use the approximation 
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mentioned above for a very large M, we can write the total Hamiltonian as 
H = He + W^'''\ We will let 

Xo = ll){X)if {xi,X2, . . .) 

be the initial state. We are interested in its time evolution and since the 
commutator [He, W^'^^'^] of the two terms appearing in the total Hamiltonian 
commutes with both, we can write [|T9|] : 

Xt = e-'*^=/^-'*^^S"'/^e'""'=^*'/^V(X)(^(xi,a;2, • • .) . (2) 

We now replace W^^^^ with the coupling W^"^^^ = 7-PI]xj, where P is the 
momentum operator of B and 7 is a constant. By following a similar path 
it is easy to write Xt on momentum basis of the macroscopic body: 

Xt = v;(P)e^""^^*'/^'e-**^^/^e-**^^S"'/V(a;i,a:2, • • • ,x„) (3) 

where ip{P) is the Fourier transform of tpiX). In the case of coupling W^^'^'^ 
we will also consider the generic entangled initial state xi, X2, ■ ■ ■ , Xn)- 
In this case: 

_ ^inakXtyn^^-itH,/h^-itkXj2x,/n^j^^^^^^^^ . . . , x„) . (4) 

We will use these equations later. 

3. THE REDUCED DENSITY MATRIX 

We will first consider the case of coupling W^^^^ and calculate the macro- 
scopic body's reduced density matrix elements Px^x" on the position basis. 
They are obtained by taking a partial trace of the total density matrix, i. 
e. integrating < X'\xt >< Xt\X" > over all the degrees of freedom of the 
environment: 

J rfa;i---da;„<^*(a;i,...,x„)^(xi,...,a;Je^*'=(^'-^")S-»M . (5) 
We used Eq. (0) and the fact that operator A = e''^^"^^ is unitary. 
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The above equation can be further manipulated by introducing, in heu of 
Xi the co-ordinate r] = J2xi/n of the centre of mass of the particle system con- 
stituting the environment and the co-ordinates of the particles with respect 
to their centre of mass: Xi = rj + ^i. The are not independent as they must 
satisfy the relation I^iLi ~ 0- ^^^^ consider 77 and the first n — 1 relative 
co-ordinates (i.e. ,^1, .^21 ■ ■ ■ j ■Cn-i) as independent variables. The last relative 
co-ordinate is expressed as ^„ = — J27=i The Jacobian determinant of the 
transformation is equal to n. We will write the integration volume element 
for the new variables as dV = drjdS, where dS — nd^id^2 ■ • • d^n-i- By 
integrating (p*{r] + ^i, . . . ,)(p{ri + $^1, . . . ,) as expressed in the new variables 
on dS and keeping r] fixed, we obtain the quantity 

w{7i) = ldScp*{rj + ^,,...,Mr^ + ^,,...,) (6) 

which is the probability density distribution of 77 in the initial state. If we 
define z — nk{X' — X")t/h we obtain 

pP^„ = / dTje'^Mv) ■ (7) 

The Fourier transform f{z) — FT[w{ri)] — J dr}e^^'^w{r]) therefore appears 
as a factor in the expression of Px^xn- In passing we note that at each 
successive instant, as is apparent from the previous expression, the absolute 
values of the off diagonal matrix elements can never be bigger than their 
initial values. If we consider the mean value rj of 77, and its mean square 
deviation a — {rj — 7})'^ (which we take as finite), for short periods of time we 
obtain: 

\f{z)\^{l-azy2) . (8) 

Now, as ^(r;) is a positive definite quantity with mean square deviation a > 0, 
in the initial moments the absolute value of the off diagonal matrix elements 
decreases, apart from the exceptional case in which a = 0. But the really 
interesting point is the asymptotic time behaviour. Here we see not only that 
^(77) is a real definite positive function, but also that its integral extended 
from minus to plus infinity must be one. So when we have an ordinary func- 
tion, it is absolutely integrable, with the result that \ f{z)\ ^ for 2; ^ ±00. 
This brings us to the important conclusion that starting from any factorized 
initial state, the off diagonal matrix elements Px'/x" (excluding exceptional 
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cases) go asymptotically to zero in the two time directions. To this end we 
should remember that if 10(7]) is m times continuously differentiable and its 
m derivatives are integrable, we obtain \f{z) \ for z ±00 more quickly 
than 1/z"^. If w{ri) is one of the more commonly encountered probability 
density distributions (such as a Lorentz or a Gaussian distribution), the off 
diagonal matrix elements tend to zero like the related Fourier transform, an 
exponential or a Gaussian, respectively. When the environment initial state 
vector is factorized in the states of the individual components, the proba- 
bility distribution w{ri) is a Gaussian under very general and physically 
acceptable conditions. 

Additional insights into the result obtained above are provided by the fol- 
lowing examples, which consider the unfavourable case of a non differentiable 
function such as 



and the case in which w{r]) is a non-ordinary and somewhat singular "func- 
tion" such as a S{ri — fj). In the first case we have 



f{z) = sm[zL]/zL = sm[nk{X' - X")Lt/h]/[nk{X' - X")Lt/h] (9) 



which allows us to define a decoherence time scale r*^**'^^ = h/[nk{X' —X")di.rj\^ 
where A 17 is the width of w{ri). In the second we have 



As we can see, the absolute value Px'x" ^'^^ ^ constant of motion but in 
spite of this the phase oscillations, which get increasingly faster, might help 
towards practically cancelling out the spatial correlations. So far we have 
assumed that environment wave function (p{xi, X2, ■ ■ ■ , Xn) is normalizable in 
the usual sense (especially with respect to the variable rj), i.e. as a square 
integrable wave function. But we also have to consider the case in which, for 
instance, the state of 1] is that of a pure plane wave. In these circumstances we 
could resort to box normalization and obtain our result simply by taking the 
limit for L going to infinity in Eq. (P). In this case decoherence is immediate. 

We will now consider the case of linear coupling W^"^^^ rather than cou- 
pling W^'^^\ Using Eq (^) and following, mutatis mutanda, the steps we took 
before, we now see that the off diagonal elements p^pl of the macroscopic 
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body's reduced density matrix on the momentum basis depends on time as: 



where y = n'y{P' — P")t/h. The value of \p^"''''^^\p'p" (as is the case for W^^^^ 
interaction) depends on the Fourier transform f{y) = FT[w{r])], which means 
that time behaviour is the same as in the previous case, replacing z with y. 

The results we obtained for both couplings have a very plausible inter- 
pretation. Looking at Eq (Q) and Eq (^ we see that the initial states of the 
environment, coupled respectively to X' and X" (Eq (^) or to P' and P" 
(Eq (|^)), simply " move apart" in Hilbert space (actually in both directions 
of time). Indeed 

e'*'^'S"'/^¥p(xi,X2,...) 

gives a momentum translation amounting to tkX' for each particle of the 
environment, i.e. to tkX'n for their centre of mass (cm.). The two wave 
packets of the cm. (coupled to X' and X" or to P' and P" respectively) 
then become separate on a momentum basis by the quantity tk\X' — X"\n 
(or t'~f\P' — P"\n) and at some time no more overlap. If the mean square 
deviation of the environment's cm. momentum is denoted as Att, we may 
suppose that decoherence is established after a time t^'^^^ = A7r/kn\X' — X"\ 
(or T^'^'^^ = An /'yn\P' — P"\ respectively). We may use the indeterminacy 
relation to write t^'^^^ = h/kn\X' — X"\Ar], thus recovering the definition 
already given of a "decoherence time" scale. 

Let's go back to the spatial coupling W^^'^^ and take a generic wave func- 
tion ^{X,xi,X2, ■ ■ ■ ,Xn) as the initial state. We assume it is square in- 
tegrable and continuous with continuous first derivatives for all variables. 
Using Eq (Hp we obtain: 



where g{z) = FT[wx'X"{v)]j with an obvious extension of the notation, is 
the Fourier transform of: 

/n— 1 n— 1 

$(X', r]Hi, • • • , CiWiX", r/+ei, • • • , r^-^ ^,)dS (12) 

1 1 

The co-ordinates transformation thus introduced does not affect the assumed 
square integrability or the analytical properties of the initial state wave func- 
tion expressed in the new variables. Evidently wx'X'iv) = \'Wx'X'{v)\ 
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wx"X"{v) = \u!x"X"{v)\ absolutely integrable in r], and therefore their 
square root will be square integrable. If we apply Schwarz inequality to the 
integration in dS, we obtain 

\wx'X"{v)\ <w¥x'{v)wxnx"iv) ■ (13) 

The right hand quantity is integrable in rj, as is clear if we again apply 
Schwarz inequality and take the properties of the two functions involved 
into account. The limitation in Eq. ([T3|) implies that wx'X"{v) is absolutely 
integrable in rj and, for the properties of the Fourier transform already used, 
the off diagonal matrix elements starting from the initial entangled state 
Px'ilc" tend to zero for z —>■ ±00, i.e. in the two time directions. 
Lastly we will consider the harmonic coupling 

H^C^^) = -i/2k ^(X - = -l/2knX^ + kXY,Xi- 1/2 ^ x\ . 

i 

The last term is absorbed in V{xi, X2, . . •), the first only contributes a phase 
factor to the matrix elements Px^x" ^^^t is the linear coupling we 

have already seen. Starting from the same initial conditions as in the case of 
W^^^\ we obtain 

(he) _ (sc) -jfen(X'2-X"2)t/[2n] 
PX'X" — PX'X"^ 

For the sake of completeness we also consider the analogous quadratic cou- 
pling : 

^^(mhc) ^ -i/2iiY,{P -uxif = -l/{2iJi)nP^ + -iPY,x,--iv/2Y,x^i . 

i 

with 7 = z///i. Here too the last term is absorbed in the potential part of the 
environment hamiltonian, the first only contributes a phase factor and the 
rest is just the coupling W'^'^'^\ Starting from the initial (factorized) state 
considered in that case: 

(mhc) _ (mc) -m(P'2_p"2)t/[2^a] /-. c\ 

fJpip/i — fJpipii^ • V / 
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(14) 



The most significant result obtained using this environment model, in 
the limit of infinite mass for the macroscopic body, has been to show that 
decoherence occurs on the position basis in the case of linear coupling W^^'^^ 
and harmonic coupling under very general initial conditions. We have only 
proven that decoherence occurs as time is going towards infinity, and we 
cannot therefore exclude the hypothesis that the spatial correlations of body 
B may increase during finite periods. However, when the initial state was 
factorized we saw that the absolute values of the off diagonal elements of 
the reduced density matrix not only go asymptotically to zero but also can 
never exceed their initial value. That particular initial condition so look as 
preferred with respect this point. In the case of this type of initial condition 
we have also showed the time dependency with which the matrix elements in 
question go asymptotically to zero. The speed at which they do so depends 
on the initial state of the environment. 

In the case of coupling W^'^'^^ we defined a time scale t^^^'^ = h/[nk{X' — 
X")A7]], where A r/ is the initial distribution width of the environment's 
centre of mass position. We may assume that decoherence is established 
when \t\ ^ T^'^^\ We have not attempted a quantitative estimate of r^'^^^ 
as our model is far from realistic. Qualitatively, however, we observed that 
as well as depending on the environment through the initial mean square 
deviation of r], it is inversely proportional to the difference X' — X", to the 
strength of coupling constant and to the number of particles n with which the 
macroscopic body interacts. This sort of dependence is certainly reasonable 
and in passing we note that decoherence occurs, of course after a greater time 
interval, even with very small coupling constants and through interaction 
with a single particle. It follows that the influence of the environment on the 
quantum state of a macroscopic body cannot be ignored and has dramatic 
effects, however weak. 

Still in the case of a factorized initial state, we saw how the preferred 
basis (i.e. on which the macroscopic body's reduced density matrix becomes 
diagonal) depends on the type of interaction with the environment. It is 
the position basis in the case of linear W^^'^'^ and harmonic coupling and the 
momentum basis in the case of W^"^'^'' coupling. The preferred basis in our 
case is the one of the operator commuting with the interaction Hamiltonian. 
The importance of an interaction commuting with the preferred basis was 



first stressed by Zurek[^. Our results so provide further confirmation of 



what is generally taken for granted, i.e. that the natural preferred basis for 
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a macroscopic body is the position basis, because the system interacts with 
the rest of the world through spatial co-ordinates. 

Another point to note is that our results are independent of whether 
or not there is an interaction potential between the particles making up the 
environment. This means, in our example at least, that decoherence need not 
be thought of as linked to a thermalization process. On the other hand, this 
could also be inferred from the case of decoherence due to scattering^, ^ Q 
and, in fact, our model describe a scattering by independent particles when 
there is no interaction between them. 

To conclude we note that in the present model, which assume the whole 
system of macroscopic body plus the environment as a closed system, deco- 
herence is not an irreversible process p2|. This is made clear in our study by 



following the evolution of the pure state of the whole system, and because the 
results obtained are asymptotically valid in the two time directions. The situ- 
ation would be quite different for a more realistic permanent flux of incoming 
particles without initial correlation. But our model allows for processes with 
suitable initial conditions (though difficult to realize in practice and which 
in our case cannot be factorized) for which the spatial correlations of the 
macroscopic body increase over a finite time. None of this prevents deco- 
herence from occurring asymptotically. In my opinion this situation calls to 
mind the spread of the wave packet of a free particle. Depending on the 
initial conditions it may even contract at the outset, but it will always end 
up expanding asymptotically. 
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